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Abstract 

The singularity structure and the corresponding asymptotic behavior of a 3-brane 
coupled to a scalar field or to a perfect fluid in a five-dimensional bulk is analyzed in 
full generality using the method of asymptotic splittings. In the case of the scalar 
field, it is shown that the collapse singularity at a finite distance from the brane 
can be avoided only at the expense of making the brane world-volume positively 
or negatively curved. In the case where the bulk field content is parametrized by 
an analogue of perfect fluid with an arbitrary equation of state P = *yp between 
the 'pressure' P and the 'density' p, our results depend crucially on the constant 
fluid parameter 7: (i) For 7 > —1/2, the flat brane solution suffers from a collapse 
singularity at finite distance, that disappears in the curved case, (ii) For 7 < — 1, 
the singularity cannot be avoided and it becomes of the big rip type for a flat brane. 
(hi) For —1 < 7 < —1/2, the surprising result is found that while the curved brane 
solution is singular, the flat brane is not, opening the possibility for a revival of 
the self-tuning proposal. 
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1 Introduction 



Some time ago, an interesting idea to address the cosmological constant problem was 
proposed, based on the so-called self-tuning mechanism [TJ |2]. The simplest model con- 
sists of a 3-brane embedded in a five- dimensional bulk, in the presence of a scalar field. 
The latter is coupled to the brane in a particular way, motivated by string theory, that 
allows a flat brane world-volume solution independently of the brane tension value. It 
was, however, realized that a singularity appears in the bulk at some finite distance from 
the brane, which can also be thought as a reservoir through which the vacuum energy 
decays. 

An obvious question is then whether the development of such a singularity is a generic 
feature of these models, or under what conditions may be avoided. In this paper, we 
investigate this question in a generalized class of models. Since in this case a general 
solution cannot be found analytically, we use a powerful tool developed a few years 
ago, called the method of asymptotic splittings, that allows one to compute all possible 
asymptotic behaviors of solutions to the equations of motion around the assumed location 
of a singularity [3] . 

As a first step in our analysis, we consider an extended version of the simplest model 
allowing for curved brane world- volume. We then show that the emergence of the finite- 
distance singularity is the only possible asymptotic behavior for a flat brane, whereas 
for a curved brane the singularity is shifted to an infinite distance. We also provide a 
detailed study of the asymptotics of this model using the method of asymptotic splittings 
expounded in [3]. A preliminary version of these results was published in [3]. 

Next, we extend the previous analysis to the case where the bulk matter is described 
by an analogue of a perfect fluid with an arbitrary equation of state P = jp. In fact, 
the case of a massless bulk scalar is a particular case of such a fluid, corresponding to 
the value 7 = 1. Here, on the one hand, we are interested in the general dynamics 
of 'evolution' of such a brane-world analogue to cosmology for arbitrary 7, in order to 
reveal the various types of singularity that may develop within a finite distance from the 
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original position of the brane, and on the other hand we seek to determine conditions 
that may lead to the avoidance of the singularities shifting them to infinite distance away 
from the brane. 

In particular, we shall show that the existence of a perfect fluid in the bulk enhances 
the dynamical possibilities of brane evolution in the fluid bulk. Such possibilities stem 
from the different possible behaviors of the fluid density and the derivative of the warp 
factor with respect to the extra dimension. The result depends crucially on the values 
of the parameter 7. 

In general, we find three regions of 7 leading to qualitatively different results: 

• In the region 7 > —1/2, the situation is very similar to the case of a massless bulk 
scalar field. Here, the main result we establish is twofold: 

- The existence of a singularity at a finite distance is unavoidable in all solutions 
with a flat brane. This confirms and extends the results of earlier works that made 
similar investigations in different models, using other methods [5j E]. 

- The singularity can be avoided (e.g. moved at infinite distance) when the 
brane becomes curved, either positively or negatively. Thus, requiring absence of 
singularity brings back the cosmological constant problem, since the brane curva- 
ture depends on its tension that receives quartically divergent quantum corrections. 

• The situation changes drastically in the region — 1 < 7 < —1/2. The curved brane 
solution becomes singular while the flat brane is regular. Thus, this region seems 
to avoid the main obstruction of the self-tuning proposal: any value of the brane 
tension is absorbed in the solution and the brane remains flat. The main question 
is then whether there is a field theory realization of such a fluid producing naturally 
an effective equation of state of this type. 

• Finally, in the region 7 < — 1, corresponding to the analog of a phantom equation 
of state, we show that it is possible for the brane to be ripped apart in as much 
the same way as in a big rip singularity. This happens only in the flat case, while 
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curved brane solutions develop 'standard' collapse singularities. No regular solution 
is found in this region. 

Besides the above regions, the values 7 = —1/2, —1 are of special significance: when 7 = 
— 1/2, we again find a regular flat brane solution with the so-called sudden behavior [7], 
as well as a non-singular curved brane, while for 7 = — 1 there is only singular curved 
solution. 

As mentioned above, it would be very interesting to understand whether there are 
field theory representations reproducing the 'exotic' regions of 7 < —1/2. Obviously the 
analogy of the perfect fluid concerning the positivity energy conditions does not seem to 
apply in this case where time is replaced by an additional space coordinate. However, 
some restrictions may be applied from usual field theory axioms. Also, the formation of 
singularities discussed here is better understood in a dynamical rather than the usual 
geometric sense met in general relativity. In the latter case, cosmological singularities 
are forming together with conjugate (or focal) points in spacetime, and for this it is 
necessary that there exists at least one timelike dimension (and any number of spacelike 
ones, greater than two). The timelike dimension forces then the geodesies to focus along 
it rather than along any of the spacelike dimensions. In the problem discussed in this 
paper, the timelike dimension is on the brane while the singularities are forming along 
spacelike dimensions in the bulk. As we show below these singularities are real in the 
dynamical sense that some component of the solution vector (a, a', p) diverges there. 
Therefore, we abandon the usual interpretation according to which the universe comes 
to an end in a finite time possibly through geodesic refocusing, and instead we study 
how dynamical effects guide our brane systems to extreme behaviors. 

The structure of this paper is as follows: In Section 2, we analyze the case of a bulk 
scalar field. We first choose appropriate variables and rewrite the basic field equations in 
the form of a dynamical system; secondly, we introduce some convenient terminology for 
the different types of singularity to be met later in our analysis. Then, in the Subsections 
2.1 and 2.2, we give the asymptotics of the models consisting of flat and curved brane 
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respectively. In Section 3, we study the case in which the bulk is filled with an analogue 
of a perfect fluid. We first derive the form of the dynamical system and single out the 
possible dominant balances, organizing centers of all the evolutionary behaviors that fully 
characterize this case. In Subsections 3.1 and 3.2, we study carefully the asymptotics 
around collapse singularities of two types, that we call type I and type II singularities, 
respectively. In Subsection 3.3, we explore the dynamics as the brane approaches a big 
rip singularity, while in Subsection 3.4 we look at a milder singularity that resembles 
in many ways the so-called sudden (non-singular) behavior introduced in Ref. |7j. In 
Subsection 3.5, we analyze the possibility of avoiding finite-distance singularities leading 
to the existence of regular brane evolution in the bulk. Finally, in Section 4 we conclude 
and also comment on possible future work in various directions, considering for instance 
other forms of matter in the bulk. In an Appendix, we briefly outline the basic steps of 
the method of asymptotic splittings. 



2 Dynamics of scalar field-brane configuration 

In this Section, we study dynamical aspects of a braneworld model (V4 = K x ^3,^4) 
consisting of a 3-brane (.M3, #3) embedded in a five-dimensional bulk space (V5, (75) with 
a scalar field minimally coupled to the bulk. The total action S to tai splits in two parts, 
namely, the bulk action Sbuik and the brane action Sbrane, 

Stotai = Sbulk + Sbrane: (2-1) 

with 

Sbuik = I (Jjp -\iy<P)^d^ (2.2) 

Sbrane = ~ J /(0)cfyt g4 , &X Y = Y*, (2.3) 

where the measures d\i gh = d 4 xdY Vdetgs, and dfi 94 = d 4 xy^detg4, Y denotes the fifth 
bulk dimension, Y* is the assumed initial position of the brane, A is a parameter defining 
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the type of scalar field 0, k\ = M^ 3 , M§ being the five- dimensional Planck mass, and 
f(4>) denotes the tension of the brane as a function of the scalar field. 

Varying the total action ( 12. ip with respect to g$, we find the five- dimensional Einstein 
field equations in the form jl], 

1 , D _x..2^ .v7. !„ ,^2\ 24 5(y/Mf 4 f ((/>)) , axfi 



6 ^ 9 (2-4) 

while the scalar field equation is obtained by variation of the action ( 12. ip with respect 
to cf) |4] and it is: 

A°50= 75== ji o{Y), (2.5) 

Vdet# 5 d(p 

where A,B = 1,2,3,4,5 and a, (3 = 1,2,3,4 while S(Y) — 1 at Y — Y* and vanishing 
everywhere else, and 



□ 5 = ^^V A ( V / d^/ B V B 0). (2.6) 
Vdet# 5 

In the following, we assume a bulk metric of the form 

g 5 = a 2 (Y)g 4 + dY 2 , (2.7) 

where g 4 is the four-dimensional flat, de Sitter or anti de Sitter metric, i.e., 

9i = -dt 2 + f k g 3 , (2.8) 

with 

g 3 = dr 2 + h 2 k g 2 , (2.9) 

and 

g 2 = d9 2 + sin 2 6dip 2 . (2.10) 

Here fk = l,cosh(Ht)/H,cos(Ht)/H (H -1 is the de Sitter curvature radius) and hk = 
r, sin r, sinh r, respectively. 
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The field equations (12.41) and (12.51) then take the form 

a \ K ^rh' 2 kH 2 

= ^- + ^-, (2.11) 
a 2 12 a 2 

- = -- (2.12) 

0" + 4-0' = 0, (2.13) 
a 

where the prime (') denotes differentiation with respect to Y, and = 0, ±1. The 
variables to be determined are a, a' and 0'. These three equations are not independent 
since Eq. ( I2.12p was derived after substitution of Eq. ( 12. lip in the field equation G aa = 



K§T aOL} a — 1,2, 3, 4, 



a' 2 kH 2 . 2 &' 2 



- + ___ = _A/c f-. (2.14) 
a a cr o 

In our analysis below we use the independent equations (12 . 12[) and (I2.13P to determine 
the unknown variables a, a' and (J)', while Eq. (12. lip will then play the role of a constraint 
equation for our system. 

Assuming a Y — > —Y symmetry and solving the Eqs. ( 12 .4p (the -aa- component, 
a = 1, 2, 3, 4) and (I2.5P on the brane we get 

a'(n) = ~f^ Y *))<y*), (2.15) 

<P\Y*) = 2A ' ^ > 

The particular coupling used in [1] allows only for flat solutions to exist. This easily 
follows by using equations (I2.15P and ( I2.16P and solving the FRW equation ( 12. lip on the 
brane for kH 2 : 

Clearly, k is identically zero if and only if: 

7W = 2 V i 145 ' 

or equivalently, if and only if /(0) oc e 2y ^^ Kr ^ (the authors of pQ have set A = 3 and 
hence the appropriate choice for the brane tension in that case is f{<p) oc e 2ft5( ^). In our 
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more general problem, the coupling function cannot be fixed this way. By working with 
other couplings we can allow for non-flat, maximally symmetric solutions to exist and 
avoid having the singularity at a finite distance away from the position of the brane. 

Our purpose is to find all possible asymptotic behaviors around the assumed position 
of a singularity, denoted by Y s , emerging from general or particular solutions of the 
system fl2.lip - fl2.13p . The most useful tool for this analysis is the method of asymptotic 
splittings [3] (see the Appendix for a brief introduction), in which we start by setting 

x — a, y = a' , z — (ft' . (2-17) 

The field equations (12.121) and (I2.13P become the following system of ordinary differential 
equations: 

x = y (2.18) 
y = -\Az 2 x (2.19) 



-Ay-, (2.20) 

x 



where A = nf/A. Hence, we have a dynamical system determined by the non-polynomial 
vector field 

f = (y,-XAz 2 x,-4y-Y . (2.21) 



x< 

Equation (12. lip does not include any terms containing derivatives with respect to Y; it 
is a constraint equation which in terms of the new variables takes the form 

y 2 AX 2 kH 2 . , 

y - 2 = — z 2 + —. 2.22 

x z 3 x z 

Equations fl2.18p - fl2.20l) and (I2.22p constitute the basic dynamical system of our study 
in this Section. 

Before we proceed with the analysis of the above system, we introduce the following 
terminology, which we use in subsequent paragraphs, for the possible singularities to 
occur at a finite-distance from the brane. Specifically, we call a state where: 

i) a — > 0, a' — > oo and p — > oo: a singularity of collapse type I. 
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ii) a — > 0, a' — > a' s and p — > 0: a singularity of collapse type Ha, 
a — » 0, a' — > a' s and p — > p s : a singularity of collapse type lib, 
a — t- 0, a' — )• a' s and p — > oo: a singularity of collapse type lie, 
where a' s , p s are non- vanishing finite constants. 

iii) a — )• oo, a' — )■ — oo and p — )■ oo: a big rip singularity, 

where p denotes the density of the matter component in the bulk that is considered 
each time. For the case of interest in this Section, p = \<f)' 2 /2 and we shall show that 
there are two major cases to be treated, the first is when we choose k = in f 1 2 . 2 2 1) 
and corresponds to a brane being flat, while in the second case k ^ 0, giving constant 
curvature to the brane. We shall treat these two cases independently. One important 
result of our analysis of this system will be that the inclusion of nonzero curvature for 
the brane moves the singularity (that is of the collapse type I class) an infinite distance 
away from the brane. 



2.1 Collapse type I singularity: Flat brane 

In this Subsection we take k = in the constraint equation (12.221) . 

V - 2 = ^ (2.23) 

We shall show that the only possible asymptotic behavior of the solutions of this system 
(flat brane) is that a — y 0, a' — y oo and <ft' — > oo, as Y — > Y s . 

Following the basic steps of the method of asymptotic splittings expounded in the 
Appendix, we start our asymptotic analysis by inserting in the system (12. 18[) - fl2.20j) the 
forms 

(x,^) = («^,/3TfVT m ), (2.24) 

where T = Y — Y s and 



(p, q, m) e Q 3 and (a, (3, () G C 3 \ {0}. (2.25) 



S 



We find that in the neighborhood of the singularity the only possible dominant balance, 
that is pairs of the form 

£? = {a, p}, where a =(«,/?, 5), p = (p, q, r), (2.26) 

determining the dominant asymptotics as we approach the singularity, is the following: 

Bt = {(a, a/4, y/3/(AVAX)), (1/4, -3/4, -1)}. (2.27) 

Note that a second balance B2 becomes only possible when we allow for nonzero curva- 
ture, k 7^ 0, and will be analyzed in the next Subsection. There are no other acceptable 
balances, hence all the possible asymptotic behaviors for a flat and curved brane can 
be described uniquely by the balances B\ and B2 respectively. Our purpose is then to 
construct asymptotic expansions of solutions in the form of a series defined by 



x = T p (a + S^ lCj T J ' /s 



(2.28) 



where x = (2, y, z), Cj = (cji, Cj 2 , Cj-3), and s is the least common multiple of the denom- 
inators of the positive /C-exponents (cf. [3], [8]). 

First we calculate the Kowalevskaya matrix (/C- matrix in short), given by 



£ = £>f(a)-diag(p), 

where -Df(a) is the Jacobian matrix of f, which in our case reads 

/ n 1 n \ 



(2.29) 



Df(x,y,z) 



\ 






1 





A\z 2 





-2AXxz 


Ayz 


Az 


_4y 


X 2 


X 


X 



J 



to be evaluated on a. For the B\ balance we have that a = 1 
p = (1/4,-3/4, —1), thus giving 

t -i 1 

3 V^AXa 



/Ci 



4 
3 



16 

V3 



4 



(2.30) 

a, a/4, v / 3/(4v / Aa)), and 
\ 

(2.31) 







9 



The next step is to calculate the /C-exponents for this balance. These exponents are 
the eigenvalues of the /Q matrix and constitute its spectrum, spec{K\). The arbitrary- 
constants of any (particular or general) solution first appear in those terms in the series 
(I2.28P whose coefficients c& have indices k = gs, where g is a non-negative /C-exponent. 
The number of non-negative /C-exponents equals therefore the number of arbitrary con- 
stants that appear in the series expansions (12.281) . There is always the —1 exponent that 
corresponds to an arbitrary constant that is the position of the singularity, Y s . The bal- 
ance Bi corresponds thus to a general solution in our case if and only if it possesses two 
non-negative /C-exponents (the third arbitrary constant is the position of the singularity, 
Y s ). Here we find 



spec(/Ci) = {-1,0,3/2} (2.32) 

so that £>i indeed corresponds to a general solution. After substituting in the system 
(12.18p - (12.20jl the series expansions 



x 



o Cjl T j ^\ y = o c j2 r^ 2 ~V\ z = E~ c j3 T j ^-\ (2.33) 



we arrive at the following asymptotic solution around the singularity: 

x = aT 1 / 4 + ^c 32 T 7 / 4 + • • • (2.34) 
y = ^r 3 / 4 + c 32 T 3 / 4 + • • • (2.35) 



4 



4VAX 7a VAX 

The last step is to check whether for each j satisfying j /2 = g with g a positive eigen- 
value, the corresponding eigenvector v of the /Ci matrix is such that the compatibility 
conditions hold, namely, 

V T . p- = 0, (2.37) 
where Pj are polynomials in c 1; . . . , Cj_i given by 

K lCj - (j/s) Cj = Pj. (2.38) 
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Here the relation j/2 = 3/2 is valid only for j = 3 and the associated eigenvector is 

aVAX 7aVAX 



The compatibility condition 



V3 ' 4V3 ' 



v T ■ {Ki - (3/2)X 3 )c 3 = 0, 



2.39) 



(2.40) 



therefore indeed holds since 

/ 



(/d - (3/2)X 3 )c 3 = c 32 



7 
4 



\ 



3 

~16 
V3 



3 
4 



\/3~AA 



v/3 



\ 4aVAX aVAX 



3 
2 



/ 



4 
7 



4v/3 



\ 7aVAX ) 





v°y 

(2.41) 

This shows that a representation of the solution asymptotically by a Puiseux series as 
given in eqs. fl2.34p - fl2.36p is valid. Hence, we conclude that at finite distance from the 
brane, a collapse type I singularity develops, i.e., as Y — > Y s the asymptotic forms of the 
variables are: 

a -»■ 0, a' -»■ oo, </>' oo. (2.42) 

This is exactly the asymptotic behavior of the solution found previously by Arkani- 
Hamed et al in pp. Our analysis shows that this is the only possible asymptotic behavior 
for a flat brane since there exist no other dominant balances in this case. 



2.2 Behavior at infinity: Curved brane 

In this Subsection we show that the collapse type I singularity that necessarily arises in 
the case of a flat brane is avoided (or shifted at an infinite distance away from the brane) 
when we consider a curved brane instead. 
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The new asymptotics follow from the study of a second balance that results from the 
substitution of (I2.24[) in (12.18p - fl2.20p . We calculate this new balance to be 



B a = {(a,a,0),(l,0,-1)}. 



(2.43) 



It corresponds to a particular solution for a curved brane since it satisfies Eq. (I2.22p for 
k 7^ and a 2 = kH 2 (here we have to sacrifice one arbitrary constant by setting it equal 
to kH 2 ), k — ±1. The /C-matrix of B2 is 

/-ll \ 




K 2 = Df(a, a, 0) - diag(l, 0, -1) 



V 







(2.44) 



3/ 



with eigenvalues 



spec(/C 2 ) = {— 1, 0, -3}. 



(2.45) 



Thus for the balance B2 we find two distinct, negative, integer /C-exponents and an 
infinite expansion in negative powers of a particular solution (recall that we had to 
sacrifice one arbitrary constant) around the presumed singularity at Y s , with the negative 
/C-exponents signaling the positions where the arbitrary constants first appear [9]. We 
therefore expand the variables in series with descending powers of T in order to meet 
the two arbitrary constants occurring for j = —1 and j = —3, i.e., 



x 



_ y-00 
, Z — 2^ j=0 Cj3 1 



(2.46) 



Substituting these series expansions back in the system f l2.18l) -( l2~20|) and after some 
manipulation, we find the following asymptotic behavior 



x = aT + c_n 
y = a H 

Z = C_33 



sT _4 + 



(2.47) 
(2.48) 
(2.49) 
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It is easy to check the compatibility conditions for j = — 1 and j = —3. We find that 

/ 1 \ / c_n \ / \ 
(/C 2 +X3)c_ 1 = 1 = , (2.50) 

V -2 j \ j \ j 



and 



/ 2 1 \ 



/ \ 





V / 



(2.51) 



(/C 2 + 3X 3 )c_ 3 = 3 

\ / \ C 33 / 

so that the compatibility conditions are indeed satisfied. The expansions given by 
Eqs. (I2.47I) - (I2.49I) are therefore valid, and we can say that as T — > 0, or equivalently 
as S = 1/T — > oo, we have that 



a — > oo, a' — > oo, 0' — > oo. 



(2.52) 



Therefore for a curved brane we find that there can be no finite-distance singularities. 
The only possible asymptotic behavior is the one given in f 1 2 . 5 2 j) which is only valid at 
an infinite distance from the brane. 

It is interesting that apart from the balances (12.271) and ( 12 .43ft . a third balance which 
initially arises from the substitution of (I2.24p in (12. 181) - (12. 201) . namely, the form 



£ 3 = {(a,0,0),(0, 



'1)}, 



is not acceptable in this case of bulk scalar field, since it does not give the necessary —1 
/C-exponent. If this balance were acceptable it would yield a finite-distance behavior of 
the type a — > a, a' — > a' s , <p' — > <p' s , where a is the constant appearing in the balance 
B3 and a' s , <fi' s are also constants. This would be similar to the sudden behavior met in 
classical four- dimensional cosmologies with a perfect fluid [7] , in the sense that the warp 
factor, in analogy with the scale factor, its derivative and the density all remain finite. 
Below, we find that such a balance does become possible when we replace the scalar field 
with a perfect fluid. 
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3 Dynamics in a perfect fluid bulk 



In this Section, we rewrite the brane model living in a bulk, which we now consider to be 
filled with a perfect fluid, as a dynamical system in three basic variables and completely 
identify the principal modes of approach to its singularities, that is we find all the 
dominant balances of the system. In this case, bulk space is filled with a perfect fluid with 
equation of state P = 7/3, where the pressure P and the density p are functions only of the 
fifth dimension, Y. We assume again a bulk metric of the form fl2.7p - fl2.10p and an energy- 
momentum tensor of the form Tab = (p + P)uaUb — Pqab, where A,B = 1,2,3,4,5 and 
ua = (0, 0, 0, 0, 1), with the 5th coordinate corresponding to Y . 
The five-dimensional Einstein equations, 



where as in the case treated previously k = 0, ±1, and the prime (') denotes differentia- 
tion with respect to Y. The equation of conservation, 



Gab — k 5Tab> 




can be written in the following form: 




(3.2) 



(3.3) 



(3.4) 



becomes 




(3.5) 



a 



Introducing the new variables 



x = a, y = a , w = p, 



(3.6) 
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Eqs. fllL2]) and take the form 



x' = y, (3.7) 

V' = -2A { -^Mwx, (3.8) 

w ' = _4(i + 7 )^ W; (3.9) 
x 



while Eq. (13. 3 p reads 



, Aw+—, A = 4/4. (3.10) 



x 2 3 x 

Since this last equation does not contain derivatives with respect to Y, it is a velocity 
independent constraint equation for the system (13 • 7f) - f T379h . 

The next step is to apply the method of asymptotic splittings in an effort to find all 
possible asymptotic behaviors of the dynamical system (I3.7p -f l3~9|) with the constraint 
(13.1 Op . by building series expansions of the solutions around the presumed position of the 
singularity at Y s . We note that the system (I3.7p -( 13T9|) is a weight homogeneous system 
determined by the vector field 

y, -2A { -^-^-wx, -4(1 + i)tu\ ^ . (3.11) 

In order to compute all possible dominant balances that describe the principal asymp- 
totics of the system, 

(2,^) = (ar^/3TUT m ), (3.12) 

where T = Y — Y s , we look for pairs of the form (I2.26P with (I2.25p . as we did in Section 2 
for the case of the bulk scalar field. We find after some calculation the following list of 
all possible balances for our basic system (13. 71) - (13. 101) : 

y Bt = | (a, ap, , (p,p - 1, -2) j ,p = * + 7 ^ -1/2, -1, (3.13) 

,B 2 = {(a,a,0),(l,0,-2)}, 7 ^ -1/2, (3.14) 

_ 1/2 S 3 = {(«,«, 0),(l,0,r)}, (3.15) 

-1/2B4 = {(", a, 5), (1,0, -2)}, (3.16) 

_ 1/2 i3 5 = {(«,0,0),(0,-l,r)}, (3.17) 
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where -i/2-Bj = 7 =-i/2 Note that, as already mentioned in the Introduction, the first 
balance 7 £>i for 7 = 1 coincides with the balance B\ in eq. ( I2.27p . where the fluid was 
described by a massless bulk scalar field with an arbitrary coupling to the brane. 

The above balances are exact solutions of the system and they must therefore also 
satisfy the constraint equation ( 13 . 1 Qj) . This fact alters the presumed generality of the 
solution represented by each one of the balances above and determines uniquely the 
type of spatial geometry that we must consider: The balances and -1/2-65 are found 
when we set k = 0, and describe a (potentially general) solution corresponding to a flat 
brane, while the balances 7 £>2 and -1/2-63 were found when k 7^ and describe particular 
solutions of curved branes (since we already have to sacrifice the arbitrary constant a by 
imposing a 2 = kH 2 ). For the balance -1/2-64, on the other hand, k is not specified and 
hence it describes a particular solution for a curved or flat brane (particularly since we 
have to set 5 = (3/(2A))(l - kH 2 /a 2 ) to satisfy eq. (13TT0|l ). 

Each one of these balances are analyzed in detail in the following Subsections accord- 
ing to the nature of asymptotic behaviors they imply. 

3.1 Collapse type I singularity 

We shall focus in this Subsection exclusively on a study of the balance 7 6i and show 
that for certain ranges of 7 it gives the generic asymptotic behavior of a flat brane to 
a singularity of collapse type I. Our analysis implies that such behavior in the case of a 
fluid bulk can only result from a 7 £>i type of balance. 

As a first step we calculate the K- matrix, /C = Df(a) — diag(p), where -Df(a) is the 
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Jacobian matrix of f \j. We have 



/ 



\ 



Df(x,y,w) 



[1 + 2 1 )Aw 



2i)Ax 



(3.1* 



V 



4(1 + 7) 



yw 

X 2 



-4(1 + 7 ) 



w 

X 



-4(1 + 7) 



y 

X 



J 



to be evaluated on a. The balance 7 $i has a = (a, ap, 3p 2 /2A), and p = (p, p — 1, —2), 
with p = 1/(2(7 + 1)). Thus the /C-matrix for this balance is 

3 



X 



m ( OP, -^jV 2 ) - di ag(p, p - 1, -2) 



— diag 



Df I (/. 
/ 



2(1 +7)' 8,4(1 + 7 ) 2 
1 



2(1 + 7) 



l + 2 7 



1 l + 2 7 

2(1 + 7)' 2(1+7)' 

N 



-2 



1 + 27 

'4(l + 7) 2 2(1 + 7) 



--{\ + 2i)Aa 







(3.19) 



/ 



\ 4(l + 7 ) 2 ,4a 2(l + 7 )Aa 
We then calculate what the /C-exponents for this balance actually are. Recall that 
these exponents are the eigenvalues of the matrix 7 /Ci and constitute its spectrum, 
spec( 7 /Ci). The balance 7 Bi corresponds to a general solution in our case if and only if 
it possesses two non-negative /C-exponents (the third arbitrary constant is the position 
of the singularity, Y s ). Here we find 

l + 2 7 



spec( 7 /Ci} 



■1,0, 



7 



(3.20) 



The last eigenvalue is a function of the 7 parameter and it is positive when either 7 < — 1, 
or 7 > —1/2. We consider here the case 7 > — 1/2 since, as it will soon follow, this range 



1 f is the vector field resulting from the dynamical system (|3.7p - (|3.9l) and {a, p} is the balance 7 2?i. 
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of 7 is adequate for the occurrence of a collapse type I singularity. The case of 7 < — 1 
leads to a big rip singularity and will be examined in Subsection 3.3. 
Let us assume 7 = —1/4 for concreteness. Then 

_ 1/4 B 1 = {(a,2Qt/3,2/(3A)),(2/3,-l/3,-2)}, (3.21) 

spec(_ 1/4 £i) = {-1,0,2/3}. (3.22) 

Substituting in the system (I3.7p -( 13T9|) the particular value 7 = —1/4 and the forms 



x 



y 



y°° r.„p/ 3-1 / 3 



W 



=0 C j3 



YJ/3-2 



we arrive at the following asymptotic expansions: 



x 



«T 2 / 3 -^c 23 T 4 / 3 + ---, 
^T- 1 /3_^ac 23 T 1 /3 + ... ; 



w 



3A 



T- 2 + c 23 T- 4 / 3 + -- 



(3.23) 

(3.24) 
(3.25) 
(3.26) 



For this to be a valid solution we need to check whether the compatibility condition holds 
true for each j satisfying j/3 — g with g a positive eigenvalue. Here the corresponding 
relation j/3 = 2/3 is valid only for j = 2 and the compatibility condition indeed holds 
since, 

\ 



U /4 Ki - (2/3)X 3 )c 2 



I 4 

3 


1 





2 


1 


Aa 


~9 


~3 


3~ 


4 


2 


2 


\ 3Aa 


Aa 


~3 



c 23 



/ Aa \ 
2~ 

2Aa 



(3.27) 



Eqs. ( l3T24j) -f l3T26|) then imply that as T ->■ 0, 



a — > 0, a' — y 00, 



p 00. (3.28) 

This asymptotic behavior corresponds to a general solution of a flat brane that is valid 
around a collapse type I singularity. We thus regain a behavior similar to the one met 
in Subsection 2.1 for the case of a flat brane in a scalar field bulk. 
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3.2 Collapse type II singularities 

In this Subsection, we show that for a curved brane (k = ±1) the long-term (distance) 
behavior of all solutions which depend on the asymptotics near finite-distance singular- 
ities turn out to be of a very different nature. In particular, we shall show that the 
balances 7 £?2 for 7 < —1/2, -1/2-63 for r < — 2 and -1/2-64 (as we have already mentioned 
7 6 2 and -1/2-63 correspond to a curved brane whereas -1/264 corresponds to a flat or 
curved brane), imply the existence of a collapse type Ha, b or c singularity. This is in 
sharp contrast to the asymptotic behavior found for a curved brane in the presence of a 
bulk scalar field (see Subsection 2.2), wherein there are no finite-distance singularities. 
For the balance 7 £>2 we find that 



7 /C 2 = Df (a, a, 0) - diag (1, 0, -2) 



/-ll \ 

-^Aa(l + 2 7 ) 
^00 -2(l + 2 7 ) J 



(3.29) 



and hence, 

spec( 7 /C 2 ) = {-1, 0, -2(1 + 2 7 )}. (3.30) 

We note that the third arbitrary constant appears at the value j = —2(1+27), 7 < — 1/2- 
After substituting the forms, 



x 



S^oC^T^ 1 , y = E? =0 c j2 Tl, w = S~ Cj3 T^ 2 , (3.31) 



in the system f l3.7p -f l3T9|) . to proceed we may try giving different values to 7: Inserting 
the value 7 = —3/4 in the system for concreteness we meet a third arbitrary constant at 
j = 1 (spec(_3/ 4 /C 2 ) = {—1,0, 1}). We then arrive at the following asymptotic forms of 
the solution: 

An 

x = aT + — c 13 T 2 + --- , (3.32) 
o 

Aa 

y = a + ~Y c i3^ H , (3.33) 

w = c 13 T- 1 + --- , (3.34) 
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where C13 7^ oB We need to check the validity of the compatibility condition for j = 1. 
But this is trivially satisfied since 

/ -2 1 \ / Aa/6 \ ( \ 
(-3/4^2 -X 3 )ci = -1 ci 3 = . (3.35) 

V J { 1 ) \o) 

The series expansions in eqs. fl3.32p - fl3.34l) are therefore valid and we conclude that 
as T — y 0, 

a 0, a' — y a, p — > 00, q^O. (3.36) 

This is a collapse type lie singularity. It will follow from the analysis below that the 
behavior of p depends on our choice of 7 (thus giving rise to three possible subcases of 
a type II singularity). Indeed, choosing for instance 7 = — 1 (spec(_i/C 2 ) = { — 1,0,2}), 
we find that the solution is given by the forms, 



Acx 3 
x = «T + — c 23 T° + 
9 

, Aa , 

w = c 23 H , 



(3.37) 

(3.38) 
(3.39) 



.1 



where C23 7^ 0[_|. Note that the compatibility condition is satisfied here as well since 



' i/Ca - 2X 3 )c 2 



We see that as T — y 0, 



1 







\ 



-2 2Aa/3 
\ 



C23 



/ Aa/9 \ 
Aa/3 

\ 1 / 



/o\ 





(3.40) 



a — y 0, a' a, p — >• c 23 , a 7^ 0. 



(3.41) 



2 If we do not set from the beginning 7 = —3/4 but instead we let 7 be arbitrary, then in the last 

step of the calculations at the j = 1 level we find that either C13 = or 7 = —3/4. 

3 Had we let 7 be arbitrary we would have found that in the step j = 2 of the procedure either c 2 3 = 

or 7 = — 1. 
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This is a collapse type lib singularity in our terminology and is clearly different from 

dSSSJ). 

A yet different behavior is met if we choose for instance 7 = —5/4. The /C-exponents 
are given by spec(_ 5 / 4 /C2) = {—1,0,3}, and the series expansions become, 



x 

y 

w 



a + — c 33 T 3 + -- 



(3.42) 

(3.43) 
(3.44) 



C33T + 

where C33 7^ 0. These expansions are valid locally around the singularity since the 
compatibility condition holds true because 

/ -4 1 \ ( Aa/12 \ ( \ 

-3 Aa 
\ ) 



(-5/4^2 - 3X 3 )Ci = 

For T — > 0, we have that 



C32 



Aa/3 
1 



/ 





\0J 



(3.45) 



(3.46) 



a — > 0, a' a, p — > 0, a ^ 0, 

which means that this is a collapse type Ha singularity. This balance therefore leads to 
the asymptotic behavior of a particular solution describing a curved brane approaching 
a collapse type II singularity, i.e. a — > and a' — )■ a. The behavior of the density of 
the perfect fluid varies dramatically: we can have an infinite density, a constant density, 
or even no flow of 'energy' at all as we approach the finite-distance singularity into the 
extra dimension at Y s , depending on the values of the 7 parameter. 

We now turn to an analysis of the balances -1/2-63, for r < —2, and -1/2B4. The 
/C- matrix for -1/2-63 is 

/-ll \ 



-1/2 



/C 3 = Df (a, a, 0) - diag(l, 0, r) 




\ -2-r 



(3.47) 



4 Here again, had we let 7 be arbitrary we would have found that in the step j = 3 of the procedure, 
either C33 = 0, or 7 = —5/4. 
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and hence, 

spec(_ 1/2 /C 3 ) = {-1, 0, -2 - r}. (3.48) 

Taking —2 — r > 0, we have two non-negative /C-exponents. (The case —2 — r < is 
considered later, in Subsection 3.5, since it is quite different, and it does not imply the 
existence of a finite-distance singularity.) For r = — 3 as an example, we substitute the 
forms 



x 



^j = C jl L J 



y 



W 



=0 C i3 



.,Ti- 3 



and arrive at the expansions 



x = aT + ■ ■ ■ , 

V = oc H , 

w = Ci 3 T~ 2 + - 

The compatibility condition is satisfied because 

/ -2 1 \ 

(-1/2^3 - 3>3)Cl = 



(3.49) 

(3.50) 
(3.51) 
(3.52) 



(3.53) 



-I ci 3 = 

v o oo/ v 1 / v°y 

and so the expansions (I3.50p - fl3.52p are valid ones in the vicinity of the singularity. The 
general behavior of the solution is then characterized by the asymptotic forms 



a —7- 0, a' — > a, p — > oo, q^0. 



(3.54) 



The balance -1/2-63 for r < —2 implies therefore the existence of a collapse type lie 
singularity during the dynamical evolution of the curved brane living (and moving) in 
this specific perfect fluid bulk. 

The balance _i/ 2 £>4, on the other hand, is one with 

/-I 1 o\ 

-1/2^4 = Di (a, a, 5) — diag(l, 0, —2) = 







V 







25_ _25 



(3.55) 



a 



a 
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and 

spec(_ 1/2 /C 4 ) = {-1,0,0}. (3.56) 

We note that the double multiplicity of the zero eigenvalue reflects the fact that there 
were already two arbitrary constants, a and 5 in this balance (recall though that S had 
to be sacrificed in order for this balance to satisfy the constraint (I3.10p ). We can thus 
write 

x = aT + • • • , (3.57) 

y = cH , (3.58) 

w = $r~ 2 + --- } (3.59) 

so that as T — )■ 0, a collapse type lie singularity develops, i.e., 

a — > 0, a' — > a, p — > oo, a^O. (3.60) 

3.3 Big rip singularities 

In this Subsection we return to the study of the balance y B\ but focus on different 
7 values. In particular, we show that when 7 < —1, a flat brane develops a big rip 
singularity at a finite distance. This new asymptotic behavior implied by the balance 
7 Bi (when 7 < —1) is equally general to the one found in Subsection 3.1. 

For purposes of illustration, let us take 7 = —2. Then the balance -2B1 and the 
_ 2 ^i-exponents read, respectively, 

„ 2 Bt = {(«,-«/2,3/(8A)),(-l/2,-3/2,-2)}, (3.61) 
spec(_ 2 /Ci) = {-1,0,3}. (3.62) 

Substituting the value 7 = — 2 in our basic system given by eqs. (13.7j) - (13.9p . and also the 
forms 

x = £°° =0 CfiT''- 1 /*, y = £°°=o c i2 T^ 3 / 2 , w = £°° =0 c j3 r^ 2 , (3.63) 
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we expect to meet the third arbitrary constant at j = 3. Indeed we find: 



x 



aT^ 2 + Uac 33 T^ + ■ ■ ■ , 

-fr- 3 / 2 + ^ac3 3 r 3 / 2 + --- 



w 



T- 2 + c 33 T + 



8A 



C33 ^ 0. 



The compatibility condition is trivially satisfied for j = 3, since the product 
3X 3 )c 3 is identically zero: 



(3.64) 
(3.65) 
(3.66) 

-2^1 - 



/ 



' 2/C1 - 3X 3 )c 3 







\ 



3 
4 



— 2Aa 
2 



( 2 - Aa ^ 



c 3 3 



/ 



Aa 



\ 



1 





W 



(3.67) 



/ 



V 4Aa 2Aa 

The series expansions given by eqs. fl3.64p - fl3.66p are therefore valid asymptotically for 
T — > so that we end up with the asymptotic forms 



a — > 00, a —7- —00, p — >■ 00. 



(3.68) 



We therefore conclude that the balance 7 £>i leads to a general solution in which a flat 
brane develops a big rip singularity after 'traveling' for a finite distance when the bulk 
perfect fluid satisfies a phantom-like equation of state, i.e. 7 < —1. Note that using 
the analogy between the warp factor of our braneworld and the scale factor of an ex- 
panding universe, we can say that this singularity bares many similarities to the one 
studied in Refs. [10], [H], [12], since it is also characterized by all quantities a, a', p, and 
consequently P, becoming asymptotically divergent. Thus, the results in this Subsection 
indicate that a flat brane traveling in a 7 < —1 fluid bulk develops a big rip singularity. 
This implements the behavior found in Subsection 3.1 of the present paper, wherein the 
same brane moving in a 7 > - 1/2 fluid bulk 'disappears' in a big bang-type singularity. 
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3.4 Sudden behavior 



As our penultimate mode of approach to the finite-distance singularity, we examine here 
the balance -1/2-65 = {(a, 0, 0), (0, —1, r)}. This balance has 



. 1/2 /C 5 = £f(a,0,0)-diag(0,-l,r) 



/ 1 \ 

1 
\0 -r / 



(3.69) 



and 



(3.70) 



spec(_i /2 £ 5 ) = {1,0, -r}, 

so we shall have to set r = 1 in order to have the necessary —1 eigenvalue corresponding 
to the arbitrary position of the 'singularity', Y s . After substitution of the forms 



(3.71) 



we find that the solution reads 

x = a + cuT + • • • , 

y = ciH , 

w = OH . 

The compatibility condition is satisfied since 

/ -1 1 \ /l\ 



(3.72) 
(3.73) 
(3.74) 



-1/2^5 — 2^) c i 




\ / 



Cll 



1 

\0J 







(3.75) 



and we see that as T — > 0, 



a — > a, a'— > Cn, p — > 0, a^0. 



(3.76) 



This clearly indicates that the brane experiences the so-called sudden behavior (cf. [7j) 
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3.5 Behavior at infinity 

A qualitatively different picture than what we have already encountered in our analysis 
of brane singularities in a fluid bulk is attained by exploiting either the balance 1 B 1 with 
— 1 < 7 < —1/2, or the balance 7 £> 2 with 7 > —1/2, or the balance _i/ 2 <B 3 with r > —2. 
We show in this Subsection that these three balances and only these offer the possibility 
of avoiding the finite-distance singularities met before and may describe the behavior of 
our model at infinity. 

We begin with the balance 7 £>i when —1 < 7 < —1/2. Choosing for instance 
7 = —4/5, we find spec(_4/s/Ci) = { — 1,0,-3}, and hence we may expand (x,y,w) 
in descending powers in order to meet the arbitrary constants appearing at j — — 1 and 
j = -3, i.e. 

x = Zj%c jl T* +!i ' 2 , y = E7~ Cj vr'"+ 8 / 2 , W = ET« CiS T'"- 2 . (3.77) 

We find: 

x = aT 5 / 2 + c_ 11 T 3 / 2 + 3/(10«)c 2 _ 11 T 1 / 2 + c„ 3 iT- 1/2 + --- , (3.78) 
y = 5a/2T 3/2 + 3/2c_! iT 1/2 + 3/(20a)c 2 l .T^ 2 - l/2c_ 3 iT" 3/2 + • • • , (3.79) 
w = 75/(8A)Tr 2 - 15/(2A*)c_i iT~ 3 + 9/(2Ax 2 )c 2 x X T- 4 + 

+ (-15/(2Aa)c_ 31 - 9/(4Aa 3 )c 3 11 ) T -5 H . (3.80) 



The compatibility conditions at j = — 1 is satisfied since 





( -3/2 


1 









/ 


1 








-4/5^1+^3)C-l = 


15/4 


-1/2 


2Aa, 


'5 


c-11 




3/2 











\ 75/(4A*) 


-15/(2Aa) 


1 






V 


-15/(2A*) 






w 



(3.81) 



26 



But for j = 



—3 we find 



-4/5^1 + 3X 3 )c_ 3 = 

/ 1/2 f \ / 

15/4 3/2 2Aa/5 

\ 75/(4Aa) -15/(2A*) 3 j 

( \ 



C-31 A 

-l/2c_ 3 i 

\ -15/(2Aa)c_ 3 i -9/{AAa 3 )c 3 _ 11 y 



-^(lOa 2 )^.!! 
\ -27/(4Aa 3 )cl 11 / 



P-3. 



(3.82) 



An eigenvector corresponding to the eigenvalue j = —3 is v T = (— 2Aa/15, Aa/15, 1), 
and hence we have 

^ T • P-3 0, (3.83) 

unless c_n = 0. In order to satisfy the compatibility condition at j = —3 we set 
c_ n = 0. The solution fl3T78|) -( l3T80|) with c_n = reads 



x = «T 5 / 2 + c_ 31 T- 1 / 2 + ---, 

j/ = 5a/2T 3 / 2 -l/2c_ 3 iT- 3 / 2 + -- 



75/(8A)T- 2 - 15/(2Aa)c_ 3 iT- 5 + 



(3.84) 
(3.85) 
(3.86) 



and it is a particular solution containing two arbitrary constants. As S = 1/T — >■ oo, we 
conclude that 



a — > oo, a — >■ oo, p — > oo, 



(3.87) 



and we can therefore avoid the finite-distance singularity in this case. 

Next we examine the balance 7 £>2 when 7 > —1/2. For 7 = 0, we have that 
spec( /C 2 ) = { — 1,0, —2}, and hence we substitute 



x 



(3.88) 
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and find: 



x = «T + c-i i — + 
y = a + A*/3c_ 23 T~ 2 H , 

23 



it; = c_ 23 T 4 + 



(3.89) 
(3.90) 
(3.91) 



The compatibility conditions at j = — 1 and j = —2 are indeed satisfied since 

to i o \ /i\ /o\ 



( /C 2 +X 3 )c_i 



1 -2A*/3 
\0 -1 J 



C-ll 







and 



/l 1 \ 



( /C 2 + 2X 3 )c_ 2 = 



2 -2Aa/3 

\ o o o y 

As^^l/T— 7>oo, we conclude that 



C- 2 3 



/ -Aa/3 \ 

V 1 / 



o 

V o / 
(°\ 



V / 



(3.92) 



(3.93) 



a — > oo, a — )• oo, p — >■ oo, 



(3.94) 



and the finite-distance singularity in shifted at an infinite distance. 

We now move on to the balance _i/ 2 £>3, r > —2. In this case we have two negative 
/C-exponents. If we choose the value r = 0, then the spectrum is found to be 



spec(_i /2 /C 3 ) = {-1,0,-2}, 



and so inserting the forms 



we obtain 



x = Zjr o c jl r> + \ y = ZjSc j2 r1, W = ET« Ci 3T'", 

x = «T + c_n, 

y = a, 

w = c_ 23 T~ 2 H , 



(3.95) 

(3.96) 

(3.97) 
(3.98) 
(3.99) 



28 



which validates the compatibility conditions 
(-1/2^3+23)0-1 = 



at j — — 1 and j = —2 since 
/ 1 \ / 1 \ / \ 



and 



(-1/2^3 + 2X 3 )c_ 2 = 
We see that as S = 1/T — > 00, 



1 
\ -1 J 

(1 1 o\ 

2 
\ j 



C-ll 



C-2 3 





v y 



V 1 / 






v y 

/o\ 

: 

W 



a' — > a, p — >■ 00, a 7^ 0, 



(3.100) 



(3.101) 



(3.102) 



a — >■ c_i 1 

so that the balance -1/2-63 for r > — 2 also offers the possibility of escaping the finite- 
distance singularities. Hence in such cases we find a regular (singularity free) evoluti 
of the brane as it travels in the bulk filled with the type of matter considered above. 
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4 Conclusions 

In this paper we studied a braneworld consisting of a 3-brane embedded in a five- 
dimensional bulk space, filled either with a scalar field or with an analogue of perfect 
fluid, giving special emphasis in the possible formation of finite-distance singularities 
away from the brane into the bulk. 

For the case of a scalar field in the bulk, we have shown that the dynamical behavior 
of the model strongly depends on the spatial geometry of the brane, in particular whether 
it is flat or not. For a flat brane the model experiences a finite-distance singularity, that 
we call a collapse type I singularity (a — > 0, a' — > 00, (j)' — > 00, as Y — > Y s ), toward 
which all the vacuum energy decays, whereas for a curved brane the model avoids the 
singularity which is now located at an infinite distance. Note that in this case, the 
collapse type I singularity is the only possible one that can develop at a finite distance 
from the (flat) brane. 
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In the second part of this paper, we studied the dynamical 'evolution' of a braneworld 
where we replaced the scalar field with a 'perfect fluid' possessing a general equation of 
state P = 7p, characterized by the constant parameter 7. For a flat brane, we find that 
it is possible to have within finite distance from the brane a collapse type I singularity 
met previously in the case of a scalar field (where this singularity was the only type 
possible, as we already mentioned above). In the fluid case, we showed that in addition 
to that singularity which appears inevitably in all flat brane solutions with 7 > —1/2, 
there are two other new types (for a flat brane): The first one is the very distinct big 
rip singularity which occurs with a — > 00, a' — > — 00, p — > 00, only when a phantom 
type equation of state with 7 < — 1 is considered. The second one is a collapse type lie 
singularity which may be described by the behavior a — > 0, a' — > a and p — > 00. This 
is less general than the collapse type I and the big rip singularities and it arises only 
when 7 = —1/2. Besides these singular solutions, we found the surprising result of flat 
branes without finite-distance singularities in the region —1 < 7 < —1/2. Moreover, 
for 7 = —1/2, there is also a solution with sudden behavior having a and a' finite and 
vanishing density p->0 0. 

In contrast to the bulk scalar field case where all curved brane solutions were regular, 
in the case of a perfect fluid in the bulk we found also singular such solutions. The 
possible corresponding finite-distance singularities are the ones comprising the collapse 
type II class. These are singularities with a — > 0, a' — > a and p — )■ 0, p s , 00 (corresponding 
to types Ha, b and c, respectively). The interesting feature of this class of singularities 
is that it allows the 'energy' leak into the extra dimension to vary and be monitored 
each time by the 7 parameter that defines the type of fluid; they all arise in the region 
7 < —1/2. On the other hand, we showed that for a curved brane the possibility of 
avoiding the finite-distance singularities that was offered in the scalar field case is still 
valid here, but only in the region 7 > —1/2. 

For illustration, we present a summary of all different behaviors we found for flat 
and curved branes in the table below, using the notation for the various singularities 
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introduced in Section 2 after Eq. (I2.22p and the balances (I3.13l) - fl3.17p . 



equation of state 


flat brane 


curved brane 


P = 19 


type 


balance 


type 


balance 


7 > -1/2 


singular type I 


-A 


regular 


7 ^2 at oo 


7 = -1/2 


singular lie 
regular sudden 


-l/2#4 
-l/2#5 


regular 
singular lie 


-1/2-63 , r > -2 at 00 
-1/2-63, r < -2 or -1/2-64 


-1 < 7 < -1/2 


regular 


7 jBi at oo 


singular IIa,b,c 


7 #2 


7 = -l 


no solution 




singular lib 




7 < -1 


singular big rip 




singular Ha 


,B 2 



An open question is whether there exist any physical constraints on 7 analogous to 
the weak and strong energy conditions of matter in ordinary perfect fluid cosmology. 
A related question is to find possible field theory realizations of the 'exotic' regions of 
7 < —1/2, where interesting solutions with unexpected behavior were found. The most 
important issue of course is to clarify the possibility of singularity avoidance at finite 
distance in flat brane solutions. There is no reason why the non-singular behavior for 
flat branes discovered here should not persist for arbitrary values of the brane tension 
and, indeed, it is to be expected that only particular asymptotic modes of behavior, i.e., 
specific detailed forms of asymptotic solutions, would depend on such values. Thus, the 
self-tuning mechanism appears to be a property of a general (non-singular) flat brane 
solution that depends on two arbitrary constants in the region — 1<7<— 1/2 (three for 
the general solution with sudden behavior when 7 = —1/2). Similarly, as we have shown 
here, the existence of singular curved brane solutions in some regions of 7 is independent 
of the sign of the scalar curvature (as long as the latter remains nonzero for curved 
branes), but the particular way of asymptotic approach to the singularity is sensitive to 
that sign and it may therefore change with different values of the brane tension. 

It would also be interesting to further investigate whether the properties of finite- 
distance singularities (and their possible avoidance) encountered here continue to emerge 
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in more general systems, such as the case in which a scalar field coexists with a perfect 
fluid in the bulk [13]. The analysis of this more involved case that allows for fluid in- 
teractions may also shed light to the factors that control how these two bulk matter 
components compete on approach to the singularity, or even predict new types of singu- 
larities that might then become feasible, as well as possible situations where they can be 
avoided. 
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A Appendix: The method of asymptotic splittings 

We refer briefly here to the basic steps of the method of asymptotic splittings. A detailed 
analysis can be found in Ref. [3]. 

Consider a system of n first order ordinary differential equations 

x = f(x), (A.l) 

where x = (x\, . . . , x n ) 6 W l , f(x) = (/i(x), . . . , / n (x)) and ' = ^p, Y being the inde- 
pendent variable. In this paper, we refrain from calling Y a time variable and giving it 
the interpretation of time. Since we are interested in singularities located at a distance 
from the brane and into the bulk, it seems more appropriate to talk about finite-distance 
singularities and give to the Y variable a spatial interpretation. The general solution of 
the above system contains n arbitrary constants and describes all possible behaviors of 
the system starting from arbitrary initial data. Any particular solution of (lA.lj) . on the 
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other hand, contains less than n arbitrary constants and describes a possible behavior 
of the system emerging from a proper subset of initial data space. 

We say that a solution of the dynamical system (1A. 1|) exhibits a finite-distance sin- 



gularity if there exists a Y s G R and a xo 6 I" such that 

lim ||x(Y;xo)|| ->• oo, (A.2) 



Y^Ys 



where || . || is any L p norm. The purpose of singularity analysis (cf. [3], [8]) is to build 
series expansions of solutions around the presumed position of a singularity at Y s in order 
to study the different asymptotic behaviors of the solutions of the system ( lA.lj) as one 



approaches this singularity. In particular, we look for series expansions of solutions that 
take the form of a Puiseux series (any log terms absent), namely, a series of the form 

x = T p (a + Ci T/ s ) , (A.3) 

where T = Y - Y s , p G Q n , s G N. 

The method of asymptotic splittings for any system of the form (lA.lj) is realized by 
taking the following steps: 

• First, we find all the possible weight-homogeneous decompositions of the vector field 
f by splitting it into components f^': 

f = f(o) + f(i)+... + f(*) j (A.4) 
with each of these components being weight homogeneous, that is to say 

f (j) (aT p ) = rP+^-iJfC^a.) j = Q,...,k, (A.5) 



where ael" and q^> are the positive non-dominant exponents that are defined by ( 1A.7I) 
below. 

• We substitute the forms x = aT p in the system x' = f(°)(x) in order to find all 
possible dominant balances, i.e. finite sets of the form {a, p}. The order of each balance 
is defined as the number of the nonzero components of a. 
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• For each of these balances we check the validity of the following dominance condi- 
tion ' 

E* fW)(aT p ) 

lim \ '- = 0, (A.6) 

and define the non-dominant exponents qV\ j = 1, . . . , k by the requirement that 

T p-i ~ T g . (A.7) 

The balances that cannot satisfy the condition ( 1A.6I) are then discarded. 

• We compute the Kovalevskaya matrix /C defined by 

/C = Df (0) (a) -diagp, (A.8) 

where Df^(a) is the Jacobian matrix of evaluated at a. 

• We calculate the spectrum of the /C-matrix, spec()C), that is the set of its n eigen- 
values also called the IC-exponents. The arbitrary constants of any particular or general 
solution first appear in those terms in the series ( 1A.3I) whose coefficients have indices 
k = gs, where g is a non- negative /C-exponent and s is the least common multiple of the 
denominators of the set consisting of the non-dominant exponents and of the positive 
/C-exponents (cf. [3], [8]). The number of non-negative /C-exponents equals therefore the 
number of arbitrary constants that appear in the series expansions of ( 1A.3I) . There is 
always the —1 exponent that corresponds to the position of the singularity, Y s . A dom- 
inant balance corresponds thus to a general solution if it possesses n — 1 non-negative 
/C-exponents (the nth arbitrary constant is the position of the singularity, Y s ). 

• We substitute the Puiseux series: 



x 



, ; -£~ c/r^/V = l,...,n, (A.9) 



in the system (1A.1I) . 

• We find the coefficients Cj by solving the recursion relations 



Kcj - 3 - c i = p i( c i- • • • > c i-i) ( A - 10 ) 
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where Pj are polynomials that are read off from the original system. 

• We verify that for every j = gs, with g a positive /C-exponent, the following 
compatibility conditions hold: 

v T -P i = 0, (A.ll) 

where v is an eigenvector associated with the positive /C-exponent g. 

• We repeat the procedure for each possible decomposition. 

We note that if the compatibility condition above (Eq. OA.lip ) is violated at some 
eigenvalue in the spec(/C), then the original Puiseux series representation of the solution 
cannot be admitted and instead we have to use a ip -series for each one of the eigenvalues 
with this property. This is a series that includes log terms of the form 

x = T p (a+ E^jCyT^T'logT)^) , (A. 12) 

where g is the /C-exponent for which the compatibility condition is violated. The rest of 
the procedure in this case is the same as before. 
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